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Abstract 

Perhaps the most important question in particle physics today is whether the boson 
with mass near 125 GeV discovered at the Large Hadron Collider (LHC) is the Higgs 
Boson of the Standard Model. Since a particularly important property of the Standard 
Model Higgs is its role in unitarizing VLlVLl scattering, we study the ability of the 
LHC to probe this process in the case of same-sign W pair production. We find that 
the use of the Matrix Element Method increases the significance with which the Higgs 
sector can be probed in this channel. In particular, it allows one to distinguish between 
a light and heavy SM Higgs in this channel alone with a high degree of significance, 
as well as to set important limits in the parameter space of the Two Higgs Doublet 
Model and the Strongly-Interacting Light Higgs Model with less than 200 fb" 1 at the 
14 TeV LHC, thus providing crucial information about the putative Higgs boson's role 
in electroweak symmetry breaking. 



Current address. 

1 



The mechanism of electroweak symmetry breaking (EWSB) can be tested most directly 
in high energy vector boson scattering. In fact, the tree-level amplitudes for scattering of 
longitudinally polarized W^W^, Z L Z L and W^Zi, pairs involving only vector bosons grow 
unboundedly with energy until they violate the unitarity limit. This unphysical growth 
must be canceled by the dynamics of EWSB. For instance, in the Standard Model (SM), 
this is achieved by diagrams involving Higgs boson exchange. After the recent discovery of 
a Higgs-like resonance with a mass of 125-126 GeV [TJ it will be essential to explicitly test 
whether this particle, other new physics, or a combination of the two are responsible for the 
unitarization of vector boson scattering. 

In pp collisions, this question can be studied through processes of the type pp —> jjV*V£ — > 
jj V1V2, where stands for a (off-shell) W- or Z-boson. Since these are four-body pro- 
cesses, the cross sections are small: even for yfs = 14 TeV, the typical signal rates are 0(fb) 
or less. Furthermore, only a fraction of this rate is contributed by longitudinal vector bosons, 
so the analysis of such processes at the LHC is very challenging. Consequently, much effort 
has been invested in studying high energy vector boson scattering and the improvement 
of signal selection cuts [2H1]. Most of these papers focused on counting signal events or 
analyzing individual distributions (such as the VV invariant mass distribution). However, 
additional information may be gleaned from various angular correlations, which are sensitive 
to the details of the unitarity-restoring dynamics. 

The Matrix Element Method (MEM) j5j|6] is a promising approach for comprehensively 
taking into account all information from an event and thus maximizing the sensitivity. It 
provides an algorithm for computing the likelihood that a given experimental event sample 
agrees with a theoretical model, using parton-level matrix elements and a set of input pa- 
rameters. For each single event, with observed momenta p" s , the likelihood that it agrees 
with a given model and set of model parameters a is defined as 



Here and /; are the parton distribution functions for the initial-state partons k and 
I, respectively, M. kl is the theoretical matrix element, and a a is the total cross section, 
computed with the same matrix element. The momenta pj of any invisible particles, such as 
neutrinos, are integrated over the available phase space. The combined logarithmic likelihood 
for a sample of N events approximately converges to a x 2 distribution: 



where p™ are the measured momenta of the nth event. 

By searching for the maximum of the likelihood, one can discriminate between several 
models and/or determine the parameters of a given model. The method is particularly 
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powerful for the measurement of processes with low event yield and/or unreconstructible 
event kinematics due to the presence of invisible particles in the final state. Typical examples 
are top quark physics at the Tevatron [SHE], Higgs searches [HUH], an d the production of 
weakly interacting dark matter particles at the LHC [Tlj|f2]. 

This letter reports on the application of the MEM to the process pp — > jjW + W + — > 
jjt r i ,+ UiVii, where = e,/i and j denotes a light quark jet. While this channel has a 
rather small cross section, due to the restriction to only positive charge and leptonic decay 
channels of the W bosons, it has the advantage of being experimentally clean and having 
low background. As a result, its sensitivity can be competitive with other vector boson 
scattering channels [HIH]. An analysis of the W + W~ channel, which has a larger signal rate, 
with the MEM is left for a future publication. 

For the likelihood calculation and the cross-section normalization in eq. (pQ), the complete 
set of tree-level diagrams for the partonic processes qq — > q'q'W + W + — > q'q'£ + £ l+ vnv^ (q, q' = 
u, d, s, c) have been included, i. e. all diagrams with on-shell intermediate W bosons. Besides 
the contribution from logitudinal W /+ V4 /+ scattering, this set also includes the irreducible 
background from all other diagrams leading to the same final state, q'q'W + W + . It has been 
shown that the difference between considering only these diagrams and considering the full 
process qq — > q' q' £ + £' + viv^i , including off-shell and t-channel W exchange, is negligible for 
high energy WW scattering |4J. The analysis is performed at the parton- level — but note that 
the inclusion of transfer functions for jet smearing and initial-state radiation is conceptually 
straightforward [6j[7J[l3] (see also [10]), although it substantially increases the computing 
time. 

Signal events are defined through a set of preselection cuts: 

p T/ > 20 GeV, 
p T j > 30 GeV, 

\Vh > 4, 

m ej > 190 GeV 

Here pt,« and rji are the transverse momentum and pseudorapidity, respectively, of a final- 
state object i = £,j; rriij is the invariant mass of the two objects % and j; and ARij = 
a/ (rji — r/j) 2 + (0j — 4>j) 2 quantifies the separation of two objects in the plane of pseudo- 
rapidity and azimuthal angle. The first two lines in eq. describe the general detector 
acceptance and object isolation cuts. The third line identifies the typical vector boson fusion 
topology, where the two jets are expected to go in the forward and backward directions, 
respectively. The last line removes background from ti production, which can produce an 
apparent same-sign lepton signal due to the small, but non- negligible, probability for lepton 
sign misidentification in the detector [1]. The invariant mass of a lepton-jet pair originat- 
ing from top quark decay is bounded from above by the top mass m t , so the requirement 
that the invariant mass of any lepton and jet is sufficiently larger than m t eliminates that 
background, while about 25% of the signal events are retained. 

The determination of the MEM weights (i. e. the numerator in eq. (CQ)) and cross section 
normalization factors (denominator in eq. (CE])) has been carried out with a specialized private 



\Vi\ < 2.5, 

\r}j\ < 5, AR jj4j>ee > 0.4 
\rjj\ > 1, m jlj2 > 100 GeV, 
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Figure 1: Statistical discrimination between different SM Higgs masses from analyzing the 
process pp — > jjW + W + — > jj£ + using the MEM (circles) vs. the di-lepton mass 
distribution (squares). The results are based on 100 events simulated for m HjVe{ = 125 GeV 
and y/s = 14 TeV. The error bars indicate the uncertainty from the event statistics. 



co using diagrams generated with FeynArts 3.3 [H]. As a cross-check, the results have 
been cross-checked against MadGraph/MadEvent 5 [To] and Mad Weight 2.5 [12J and 
good agreement has been founcjfl. MadEvent has also been used for the generation of the 
simulated "experimental" events that are fed into the MEM fit. Throughout this paper, the 
"experimental" events sample is based on the SM with m H>rei = 125 GeV. The cuts ([3]) have 
been consistently applied both to the event generation and weight normalization. 

In the following the analysis of the process pp — > jjW + W + — > with the 

MEM is presented for three characteristic models: the Standard Model (SM), the Two 
Higgs Doublet Model (THDM), and the Strongly Interacting Light Higgs (SILH) Model. 



SM: For m# = 125 GeV, the signal cross section for a/s = 14 TeV after the preselection 
cuts is 0.590 fb, resulting in a signal yield of 100 events with an integrated luminosity of 
170 fb -1 . Assuming this number of signal events, the results of the MEM likelihood fit are 
shown in Fig. [TJ Specifically, if one wants to test whether the unitarization of WW scattering 
is facilitated by a light Higgs boson with ran — 125 GeV or a heavy Higgs boson with 
ran = 1000 GeV, these two hypotheses could be distinguished with a statistical significance 
of more than three standard deviations. 

For comparison, Fig. [1] shows the statistical discrimination obtained from analyzing the 
m« distribution as suggested e.g. in Ref. [Bj. When using two bins in the range m« £ 
[0, 1000] GeV and a sample of 100 events, the significance stays below one standard deviation 

2 The code is available upon request from the authors. 

3 It was difficult to reach the required precision for the cross section values with MadGraph/MadEvent 
within a reasonable amount of computing time, so the results presented here are based on our private code. 
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for Higgs masses up to 1 TeV. With larger numbers of bins the discriminative power is even 
lower. 

Of course, in light of electroweak precision tests and the evidence for a 125-GeV scalar 
at LHC [H], the SM with a heavy Higgs boson is effectively excluded. Nevertheless it is illus- 
trative to discuss this scenario as an simple example that tangibly highlights the difference 
between the MEM and a traditional analysis strategy. 



THDM: A more realistic scenario is given by the THDM with a light CP-even Higgs boson 
(h°) of mass rrih = 125 GeV, which explains the observed Higgs signal, and a heavy CP-even 
Higgs {H°) with unconstrained mass uih [16J. Assuming CP conservation, the two physical 
states h° and H° are mixtures of the CP-even components of the two Higgs doublets, H® 



and Hi: 



h° = cos a H® — sin a H®, 
H° = sina H° + cosaH°. 



(4) 



Denoting the ratio of the vacuum expectation values by (H®)/ (H®) = tan/3, the Higgs- W-W 
couplings read 

g(h°WW) T}lDM 

cos(p — a) = cos£, 



g(HWW) SM 
g(H°WW) T]1DM 



(5) 



sin(/3 — a) = sin£, 



g(HWW) SM 

where g(HWW)sM is the HWW coupling in the SM. If £ = (3 — a is non-negligible, both 
h° and H° play a role in the unitarization of WW scattering. As shown in Fig. |2j the 
interplay of the two Higgs bosons in this process can be tested experimentally with the help 
of the MEM. In particular, a sample of just 100 events will be sufficient to rule out the 
parameter region with large values of £ > 7r/4 and m H o > 600 GeV at about 90% confidence 
level [assuming the data agrees with the SM]. 



SILH: The SILH paradigm encompasses a class of models with strong dynamics at the 
scale A ~ 4nf > 1 TeV and a light composite Higgs boson, which is a pseudo-Goldstone 
boson of some global symmetry [IT]. As a low-energy effective theory, it contains a light 
Higgs whose coupling to W and Z are modified by a factor 1/ ^/l — cv 2 / f 2 , where c is a 0(1) 
number. As a result, unitarization of high energy vector boson scattering is not achieved 
by the light Higgs alone, but requires the presence of additional heavy scalar and vector 
resonances, which emerge from the strong sector. Here it is assumed that these resonance 
are beyond the reach of the LHC, so that the only observable effect are the modified hWW 
and hZZ couplings^, which is equivalent to the limit m H a — > oo of the THDM. 

Similar to the previous examples, the deformation parameter cv 2 /f 2 can be constrained 
by analyzing high energy W + W + through the process pp — > jjW + W + — > jji^ 'i ,Jr vgv^ . The 
output of the MEM as a function of this parameter is shown in Fig. [31 together with the 

4 The explicit inclusion of the heavy resonances in the MEM does not pose any conceptual problem, but 
due to the substantial amount of computing time involved this is left for future work. 
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Figure 2: Projected constraints on the THDM from the MEM analysis of the process pp — > 
jjW + W + — > jjl + £' + veV£>, as function of the heavy Higgs mass m H o and the mixing angle £ 
(the light Higgs mass is fixed to m/jO = 125 GeV). The contour lines indicate A% 2 = 1,4,9. 
The results are based on 100 events simulated for £ = and y/s = 14 TeV. 

results obtained from analyzing the mu distribution. It has been checked that the latter are 
compatible with the numbers in Tab. 14 of Ref. [3] within statistical errors. As evident from 
the figure, the MEM leads to an improvement of the sensitivity by more than one standard 
deviation. 

Conclusions: High energy vector boson scattering provides a unique window into the 
mechanism of electroweak symmetry breaking, but it is difficult to analyze experimentally 
at the LHC due to the small event yield. The Matrix Element Method (MEM), which is 
an automated likelihood technique incorporating all relevant event and theory information, 
significantly improves the sensitivity for this process as compared with traditional analysis 
methods. This had been demonstrated explicity here for high energy W + W + scattering. 
As concrete examples, the method has been applied to three characteristic examples, the 
Standard Model with variable Higgs mass, the Two Higgs Doublet Model, and the Strongly 
Interacting Light Higgs Model, but it can be adapted straightforwardly to other models of 
electroweak symmetry breaking by implementing the relevant matrix elements. 

It is worth pointing out that the method does not rely on the observation of a resonance. 
In fact, for high energy same-sign WW scattering, the Higgs boson or any other unitarity- 
restroring particle contributes only in the t-channel. Consequently, the MEM will be useful 
even for a so-called "nightmare" scenario with very broad resonances |18j . 

The results presented here are based on a parton-level analysis. For a more realistic 
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Figure 3: Statistical discrimination between different values of the deformation parameter 
cv 2 /f 2 in SILH models from analyzing the process pp — > jjW + W + — > jj£ + £' + viVii using the 
MEM (circles) vs. the di-lepton mass distribution (squares). The results are based on 100 
events simulated for m# ire f = 125 GeV, ch = 0, and = 14 TeV. The error bars indicate 
the uncertainty from the event statistics. 



picture, one needs to consider systematic uncertainties and detector effects. The largest 
systematic error is related to the measurement of the jet energy, which can be taken into 
account by incorporating jet smearing functions and by treating the overall jet energy scale 
as a free parameter in the fit j6]. While this leads to a substantial increase in computing 
time, the senstivity of the MEM is not significantly reduced. Other systematic effects include 
next-to-leading order QCD corrections, which are small for vector boson fusion processes [T9] . 
and uncertainties in the quark and antiquark parton distribution functions, which can be 
controlled using data on Drell-Yan production. Therefore it is expected that the effectiveness 
of the MEM will not be significantly reduced by systematic errors. 
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